A combinatorial approach to metamaterials discovery 
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We report a high through-put combinatorial approach to photonic metamaterial optimization. 
The new approach is based on parallel synthesis and consecutive optical characterization of large 
numbers of spatially addressable nano-fabricated metamaterial samples (libraries) with quasi- 
continuous variation of design parameters under real manufacturing conditions. We illustrate this 
method for Fano-resonance plasmonic nanostructures arriving at explicit recipes for high quality 
factors needed for switching and sensing applications. 



Some fifteen years ago a paper reporting a combinato- 
rial approach to materials discoveries revolutionized ma- 
terials research and other disciplines such as chemistry 
and pharmacology by demonstrating a method for high- 
throughput parallel synthesis and analysis of novel arti- 
ficial chemical compounds [l], Q . Here we report on the 
first application of the combinatorial approach to dis- 
covery of electromagnetic metamaterials and their opti- 
mization. Electromagnetic metamaterials are manmade 
media with all sorts of unusual functionalities that can be 
achieved by artificial structuring on a smaller length scale 
than that of the external stimulus. We apply the com- 
binatorial approach to the optimization of an important 
class of plasmonic metamaterials supporting Fano reso- 
nances [3] that recently became a prime platform for new 
switching, gain and sensing applications, slow light and 
polarization control devices (for a review see Q). These 
metamaterials are arrays of asymmetric split ring metal- 
lic wire resonators or complementary arrays of ring slits 
in a metal film. Their application potential depends on 
the presence and characteristics of reflection and trans- 
mission peaks associated with a particular mode of elec- 
tromagnetic excitation that is weakly coupled to free- 
space (closed or trapped mode). The characteristics of 
this resonance are very sensitive to the environment and 
are responsive to its changes which makes them a prime 
choice for active metamaterial applications. These reso- 
nances can be tailored by design: their spectral position, 
width and depth depend on many parameters, most no- 
tably on the type of metal used for their fabrication, size 
of the ring and characteristics of the split Opti- 
mization of resonance characteristics of these metama- 
terial structures is relatively straightforward in the mi- 
crowave and terahertz parts of the spectrum using nu- 
merical electromagnetic modeling tools. For the optical 
part of the spectrum when the meta-molecules should be 
structured on the nanoscale numerical modeling becomes 
considerably less reliable for two reasons. First, accurate 
reproduction of the idealized design parameters is not 
possible due to limitations of the nanofabrication tech- 
nologies, in particular sharp edges and ideal vertical cuts 
are not possible to achieve while the surface of metals 



has roughness on a scale comparable with the smallest 
features of the design. Second, electromagnetic mate- 
rial properties used in calculations, in particular those of 
nanostructured metals are not known to the necessary 
accuracy and are normally replaced by those of bulk me- 
dia. However, they can be drastically different from those 
of bulk metals due to the higher role of surface electron 
and defect scattering in granulated and nanostructured 
metals. The unsatisfactory character of this substitute 
is well recognized by the metamaterials research com- 
munity where in some calculations arbitrarily corrected 
values of bulk dielectric parameters are used for nanos- 
tructured metals to achieve a better fit to experimental 
data. 

The combinatorial approached presented below aims 
to address the difficulty of optimizing photonic metama- 
terial designs through computer modeling by resorting to 
a parallel synthesis and consecutive optical characteriza- 
tion of a large number of spatially addressable nanofabri- 
cated metamaterial samples [libraries, see Fig. Ufa)] with 
quasi-continuous variation of design parameters. We also 
compare the results of combinatorial optimization with 
those of numerical modeling. 

We studied metamaterials consisting of arrays of 
square split-ring aperture meta-molecules in the field of 
parameters defined by the overall size s of the metama- 
terial unit cell and the asymmetry of the split [see Fig. 
DIb)]. Here s x s is the overall size of the square unit 
cell while in all cases the split ring aperture's outer and 
inner sides a and I are 100 nm and 200 nm smaller than 
the unit cell, respectively. The design slit width w was 
50 nm and parameters t and b define the asymmetry of 
the split and its size. Two libraries of regular 30 x 30/im 2 
arrays of split ring apertures in a 30 nm thick layer of gold 
were fabricated by electron beam lithography on a glass 
wafer. In both libraries the unit cell size s was varied 
from 400 nm to 500 nm in steps of 25 nm. In library 
A parameter b — was fixed and a gap size asymmetry 
a = t/l ranging from 0.1 to 0.7 in steps of 0.1 was intro- 
duced. In library B the gap size was fixed at 0.3/ and a 
gap position asymmetry /3 = (t — b) /(t + b) ranging from 
to 1 in steps of 1/7 was introduced. 
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FIG. 1: (Color online) Combinatorial library, met a- molecules and resonant modes, (a) Schematic of a combinatorial 
metamaterial library, (b) Schematic of the metamaterial unit cell, (c) and (d) show scanning electron micrographs of meta- 
material unit cells of size s = 500 nm with small and large gap position asymmetries, f3 = 1/7 and 1. (e) Characteristic I 
symmetric and II-III antisymmetric modes of excitation associated with the absorption resonances for asymmetries /3 = 1/7 
and 1. The plotted field distributions correspond to excitation wavelengths of 820, 1025 and 1200 nm, respectively. 



The transmission and reflection characteristics of these 
photonic metamaterials were measured in the spectral 
range from 800 to 2000 nm for waves polarized parallel to 
t using a microspectrophotometer. These properties were 
also simulated with a full three-dimensional Maxwell fi- 
nite element method solver in the frequency domain using 
the Drude model for the dielectric constant of gold [9[ . 

Figure EJa)-(c) illustrates for a fixed meta- molecule 
size of s = 500 nm how the metamaterial's transmission, 
reflection and absorption properties depend on the gap 
size asymmetry a = t/l when 6 = 0. In the symmetric 
case the metamaterial has broad reflection and transmis- 
sion maxima around 800 nm and 1400 nm respectively, 
while the absorption spectrum is featureless. Symme- 
try breaking, a > 0, leads to the appearance of a Fano 
resonance near 900 nm which is associated with a nar- 
row peak in absorption. Furthermore, at slightly shorter 
wavelengths, a window of transparency opens up, which 
may be discussed in terms of electromagnetically induced 
transparency |10| . Larger symmetry breaking leads to 
broadening of the absorption resonance, which red-shifts 
due to the increasing size of the top aperture. 

Figure [3^a) -(c) shows how the metamaterial transmis- 
sion, reflection and absorption depend on the gap posi- 
tion asymmetry j3 = (t — b)/(t + b) for a fixed 90 nm gap 
and a unit cell size of s = 500 nm. Similarly to the gap 
size asymmetry discussed above, a gap position asymme- 
try f3 > leads to the appearance of a Fano resonance 
at 1000 nm which broadens and red-shifts with increas- 
ing asymmetry. Intriguingly, simulations reveal an ad- 
ditional narrow absorption resonance at 830 nm, which 
should be supported by the symmetric structure and van- 
ishes with increasing asymmetry f3. The fundamentally 
different nature of these resonances becomes clear when 



comparing the associated modes of excitation shown in 
Fig. Die). The symmetric structure supports a sym- 
metric mode I, where top and bottom apertures interact 
through in-phase magnetic fields. Symmetry breaking al- 
lows the excitation of an anti-symmetric mode II, where 
the apertures of each split ring interact through anti- 
phase magnetic fields. These modes can be understood as 
the hybrid modes of two coupled resonators (apertures): 
I being the higher energy anti-bonding and II being the 
lower energy bonding mode (ll) . With increasing sym- 
metry breaking the bonding mode II becomes dominated 
by the longer aperture, see III. 

As discussed in the introduction, narrow Fano reso- 
nances have a wide range of applications. To quantify 
the width of the bonding resonance we introduce a fig- 
ure of merit defined as FOM = Ao/AA, where Ao is the 
resonant wavelength and A A is the full width half maxi- 
mum of the absorption peak. As illustrated by Figs. [2fd) 
and[3Kd), for both types of symmetry breaking smaller 
asymmetries lead to narrower Fano resonances, while the 
resonance width does not strongly depend on the meta- 
molecule size, which controls the spectral position of the 
resonance. 

As illustrated by Figs. [2]and[3j numerical results and 
experiments show similar qualitative behavior, indicating 
that simulations do provide useful guidance for the design 
and understanding of photonic metamaterials. However, 
quantitative comparison shows substantial differences be- 
tween experimental and numerical results. 

In particular, the spectral position of the measured ab- 
sorption resonances is red-shifted by about 80 nm relative 
to the simulation results. This shift may be explained 
by uncertainty regarding the real part of the dielectric 
constant of nanoscale volumes of gold at optical frequen- 
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FIG. 2: (Color online) Library A - gap size asymmetry a — t/l, where parameters I = s — 200 nm and 6 = remain fixed, 
(a) Transmission, (b) reflection and (c) absorption spectra as a function of symmetry breaking a for unit cell size s = 500 nm. 
The dashed line corresponds to the absorption maximum, (d) Combinatorial maps of resonance sharpness as a function of unit 
cell size s versus gap size asymmetry a. Here the resonance width is measured as a figure of merit, FOM = Ao/AA, defined 
as the resonance wavelength Ao of the absorption peak divided by its full width half maximum A A. Note that the experiments 
(top) cover a slightly smaller range of asymmetries than the simulations (bottom). 
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FIG. 3: (Color online) Library B - gap position asymmetry /3 — (t — b)/(t + b), where parameters I = s — 200 nm and 
t + b = 0.7/ remain fixed, (a) Transmission, (b) reflection and (c) absorption spectra as a function of symmetry breaking 
f3 for unit cell size s = 500 nm. The dashed line corresponds to the absorption maximum, (d) Combinatorial maps of 
resonance sharpness FOM as a function of unit cell size s versus gap position asymmetry j3. The resonant modes of excitation 
corresponding to I-III are shown in Fig. [TJe). 



cies. Also systematic manufacturing inaccuracies like the 
inability to fabricate sharp corners at the nanoscale con- 
tribute to a spectral shift of resonances, see Fig. QJc)-(d). 
In general spectral positions of Fano-type resonances in 
plasmonic metamaterials are also extremely sensitive to 
thin cover layers, or layers of oxides or contamination, as 
has recently been demonstrated by a 3-fold increase in 
metamaterial transmission resulting from application of 
a single layer of graphene [12| . 

While simulations predict remarkably narrow reso- 
nances of up to FOM = 50 for small asymmetries around 



0.1, the experimentally observed resonances are broad- 
ened FOM < 12 and weakened so much, that they can- 
not be clearly identified for asymmetries a < 0.3 and 
j3 < 3/7. Small size variations and manufacturing tol- 
erances in rep roducing a perfect square array of meta- 
molecules [13| , rough gold surfaces as well as finite array 
sizes Q were shown to lead to resonance broadening and 
weakening in real structures, while these factors are not 
accounted for in simulations of infinite arrays of identical, 
ideally smooth meta- molecules. Additionally, the width 
of resonances strongly depends on the imaginary part of 
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the dielectric constant of gold, which may be larger than 
the commonly used literature values [9] for nanoscale vol- 
umes of gold in the optical part of the spectrum due to 
increased rates of electron scattering. These factors can 
also explain why the remarkably narrow anti-bonding res- 
onance associated with the symmetric mode of excitation 
has not been observed in real structures. Furthermore, 
the above factors can also cause the high background of 
absorption observed in our real metamaterial samples. 

For practical applications a trade-off between narrow 
Fano resonances (small asymmetry) and strong Fano res- 
onances (large asymmetry) must be considered. The op- 
timal design parameters (i.e. asymmetry) will always de- 
pend on the accuracy of the available nanostructuring 
technology and can be found through the high through- 
put screening method demonstrated here. In our case 
asymmetries around a, j3 c± 0.4 — 0.5 are the best choice 
for most sensor and switching applications, which is in 
stark contrast to simulations suggesting optimal asym- 
metries of around 0.1. 

In summary we report a high through-put combinato- 
rial approach to the design and optimization of photonic 
metamaterials. We show for Fano resonances in split ring 
aperture arrays, that the optimal metamaterial param- 
eters can be found through parallel synthesis and char- 
acterization of material libraries with quasi-continuous 
variation of design parameters under real manufacturing 
conditions. Being unable to take the inevitable short- 
comings of nanoscale manufacturing into account, sim- 
ulations can only provide qualitative insights into meta- 
material properties and underlying mechanisms. Here we 
have identified in detail how the asymmetry of split ring 
apertures controls Fano resonances in photonic metama- 
terials, which are relevant to various applications ranging 
from sensing, switches and slow-light devices to the "las- 
ing spaser" [15j. 
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